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SUMMARY 


The particular integrals of the second-order and third-order Prandtl- 
Busemann iteration equations for the flow of a compressible fluid are 
obtained by means of the method in which the complex conjugate -variables z 
and z are utilized as the independent variables of the analysis. The 
assumption is made that the Prandtl-GIauert solution of the linearized or 
first-order iteration equation for the two-dimensional flow of a com- 
pressible fluid is known. The forms of the particular integrals, derived 
for subsonic flow, are readily adapted to supersonic flows with only a 
change in sign of one of the parameters of the problem. 


INTRODUCTION 


For- the past several years iteration methods have been increasingly 
applied to the solution of compressible -flow problems. . The most useful 
one from the point of view of aeronautical applications and the one 
discussed in this paper is based on small, perturbations with respect to 
the undisturbed flow. The Prandtl-GIauert and Ackeret solutions in two- 
dimensional subsonic and supersonic flow, respectively, obtained by 
means of the linearization of the fundamental nonlinear differential 
equation for compressible flow, are presumed. to be known a nd are taken 
as the initial steps in this iteration process. Higher-order solutions 
are then obtained by retaining appropriate powers and products of the 
perturbation quantities. This method of iteration has been variously 
labeled the Ackeret iteration process and the Prandtl-Busemann small 
perturbation method when limited to two-dimensional subsonic flow. The 
procedure has been extended in recent years to both two-dimensional and 
axisymmetrical supersonic-flow problems . 

In a recent publication (reference l). Van Dyke succeeded in 
obtaining by trial the particular integral of the nonhomogeneous second- 
order iteration equation for the velocity potential in supersonic flow. 
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The general solution is then easily obtained by adding solutions of the 
homogeneous equation with proper regard to the bou n da r y conditions at 
the surface of the solid and at infinity. 

The purpose of the present paper is to show a procedure by means 
of which the particular integrals of the hi gher-than-f irst- order iteration 
equations c an be derived in a systematic manner. The explicit expression 
obtained for the particular integrals of the second- and third-order 
iteration equations are believed to yield essentially the solution of 
the problem of hi gh subsonic flow past an arbitrary two-dimensional 
profile, since it is never a difficult problem to supply the solutions 
of the h om ogeneous equation necessary for the fulfillment of the boundary 
conditions. It is noteworthy that the particular integrals, derived for 
subsonic flow, can be adapted to supersonic flow with simply a change 
in sign of one of the parameters. 


FUNDAMENTAL EQUATIONS 

The fundamental nonlinear differential equation governing the flow 
of a compressible fluid is 

( C 2 - u 2)|£ + ( C 2 - v2)|| - uv^ + “ 0 (1) 

where 

X, Y rectangular Cartesian coordinates in flow plane 

u, v fluid velocity comppnents along X- and Y-axis, respectively 

c local speed of sound 

The condition for irrotational motion is that 


du dv 

3? = 3x 

and leads to a velocity potential. <5 defined by 


u 



v = 3 y 


( 2 ) 
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If the "body is held fixed in a uniform stream of velocity U, the 
relation between the local speed of sound c and the speed of the 

fluid \/u 2 + v 2 • is given for adiabatic processes by 


5 = <3) 

where 

c ra sound speed in undisturbed fluid 

7 ratio of specific heats at constant pressure and constant volume 

M ra Mach number of undisturbed stream (U/c^) 

With the introduction of a characteristic length Z as unit of 
length and the undisturbed stream velocity U as unit of velocity, the 
quantities X, Y, u, v, and $ for the remainder of the analysis 
denote, respectively, the nondimensional quantities X/Z, Y /z, u/u, 
v/U, and <J>/UZ, while c and c^ retain their original meanings. By 
means of equation ( 2 ), equations (l) and ( 3 ) then become, respectively. 


and 
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- <&xx + 


- Me^vlo^Y 


= 0 


w 


~ = 1 + 7 ~ h i 2 ji _ ( u 2 + v2)J (5) 


where the subscripts X and Y denote partial differentiations with 
respect to the designated variables. 

In order to obtain the iteration equations based on small perturbations 
of the undisturbed stream, the assumption is. made that the velocity 
potential $ can be expanded in the form 


G> = X + + $2 + $3 + . . . 


(6) 
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For the purpose of defining and controlling the iteration procedure, 
the function ® n+ -j_ and its derivatives are then regarded as small 

compared with the preceding approximation $ no and its derivatives. 
From equations ( 2 ) and ( 6 ), 


u = = 1 + y + ® 2 X ® 3 X ^ • • • 


and 


v = $y = ®-jjy + $2Y + ® 3 Y + • • • 

When these expressions for u, v are introduced into equation ( 4 ), 
together with the expression for c 2 /^ 2 given hy equation (5), and 
the powers and products of and their derivatives are grouped 
according to the assumptions of the smal l perturbation method, the 
following iteration equations for the first three approximations 
4^, and $3 result: 


(i ' M «) 4 mt + - 0 


( 7 ) 


( 1 " Ho 2 )"^xx + *^yy ■ M «, 2 [(7 + 1 )®ix*ixx + (7 - 1 )®ix 0 iyy + ^iy^ixy] 


(8) 


** 

( 1 " M ~ 2 ) ®3XX + $ 3YY " M m 2< (l^ix 2 + $ 2x) + 1 ) $ 1XX + (7 - 1 ) <i ’iyx] + 


- 1 ) $ 1 XX + (7 + + ® 1 X [(7 + ^O^XX 

(7 - 1)®2Yy] + 2 (^1X®1Y®1XY + $ 2Y\xi + $ lY®2XY^)j 

( 9 ) 


For slender bodies, the first few steps of this iteration process may be 
expected to yield an accurate result with the exception of a small region 
in the neighborhood of a stagnation point. Even at stagnation points 
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the iteration method has "been shown to represent correctly the effect of 
compressibility (reference 2). The accuracy of the calculations 
obviously depends upon the number of terms determined, each additional 
term reducing the region of inaccuracy in the neighborhood of a stagnation 
point. 

The iteration equations (7) > (8), and (9) may be put into more 
familiar forms by the introduction of a new set of independent variables x 
and' y, where 



Thus, for Mo,, < 1, equation (7) is transformed into a Laplace equation, 
whereas equations (8) and (9) are transformed into Poisson equations 
with the right-hand sides composed of, respectively, double products and 
triple products of previously determined perturbation quantities. It 
is further assumed that the solution of equation (7) is available. This 
initial step in the approximation to the exact nonlinear solution is 
usually easily obtained, as it represents the Prandtl-Glauert approximation 
(reference 3 } appendix B). The purpose of the present paper is then to 
derive explicit expressions for the particular solutions of the second- 
and third-order iteration equations (8) and (9). 


CALCULATION OF THE PARTICULAR INTEGRAL OF THE SECOND-ORDER 

ITERATION EQUATION 


By introducing the independent variables x and y defined by 
equation (10), the second-order iteration equation (8) becomes 

®2xx + $2 yy = ^oo 2 j^ 1 + cr )^lx <I lxx + ®ly®lxyj (H) 

where 

a = 7 + 1 M °° 2 

2 e 2 

3 2 = 1 - 
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and where use has been made of Laplace's equation 


$ lxx + $ lyy = 0 


The procedure for obt aining the particular integrals of the higher- 
order iteration equations is based on the use of the complex conjugate 
variables z and z as independent variables. Thus, 


z = x + iy 

o 

2 = x - iy 

and the equivalence of operators 

a _ 9 a 

Sc “ ciz + S2 

a . i fa _ a \ 

SF AS sj 

a 2 .a 2 .oS 2 a 2 

±_ ,a 2 aL 

^ - 3z 2 +2 ss - 

Then equation (7) for $-[_ becomes 

a 2 ®! 

“SS' 0 (“) 

The most general real solution of this equation is 

\ = + 


( 13 ) 
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or 


= R.P.w-j_(z) = R.P.w 1 (z) 


where w^(z) is an arbitrary analytic function of z, 'Bi(z) is its 
conjugate complex, and where the symbol R.P. stands for "real port of". 
The imaginary part of wi(z) is a function say, related to 
by means of the Cauchy-Riemann equations and hence also satisfies Laplace’s 
equation. The function >}q does not represent the stream function of 

the actual compressible flow and does not appear in the final expressions 
of the particular integrals. The following relations will be found use- 
ful and are easily verified: 


$ 1 

- i(“l + 

w-L) 

*lx 

= |(w lz + W lz ) 

®iy 

= |( w lz ■ 

- *lz) 

^lxy 

=. 2^ w lzz 

“ ^lzz) 

*lxx 

11 

IOJH 

H* 

IS) 

N 

+ ^lzz) 


Then 


23 lx <I lxx = R.P- (viz + vi£)wizz 


2< *Ty^lxy = _ R.P* (viz - wiz) w lzz 


and equation (11) for the second approx ima tion $2 becomes 
°2zz = TjM m 2 R.P.[owi z Wi zz + (2 + aj^-w^z] 

= R - p -[|( w lz 2 ) z + (2 + a)(wi z wi5) z ] 
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If <&2 is defined to be the real part of a nonanalytic function vq (z,z), 
then 


w 2zz = 



(2 + a)(w lz w l5 jj 


(14) 


This equation can he integrated immediately by inspection and yields 
the general solution 


*2 = ^[f^lz 2 + ( 2 + ff )*l w lz + F ( 2 )] ( 1 5) 

where, because only the real part $2 °i w 2 i s interest, only one 

arbitrary analy tic function F(z) need be Included. The function F(z) 
satisfies Laplace's equation and is so chosen as to satisfy the required 
boundary conditions at the surface of the body and at infinity. The 
part of the expression on the right-hand side of equation (15) * excluding 
the arbitrary function F(z), is the particular integral of equation (14) 
and may be expressed in real form in the following manner: 

Suppose 


F(z) = - |zw lz 2 + (j2 + ctJw-lW^ + f(z) (16) 


where f(z) is again an arbitrary analytic function of z. 
Then with the aid of the relation 

w iz = " *>* + = $lx " i$1 y 

where use has been made of the Cauchy-Riemann conditions 


®lx = *ly 


ly 
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the expression for $2 > obtained from equation (15), becomes 

$2 = [- ^lx^ly + 2 ( 2 + a )*l $ lx + R-P-f(z)] (IT) 

The expression on the right-hand side of this equation, excluding 
R.P. f(z), namely. 


® 2 = 


(l + f)®!' - *1, 


( 18 ) 


corresponds precisely to the particular integral obtained by Van Dyke 
(reference 1) for two-dimensional supersonic flow with cr replaced 
by -cr, where for supersonic flow the definition of cr is 


cr 


y + 1 

2 



It is rather noteworthy that the particular integral of the second-order 
iteration equation (8) can be obtained for both subsonic and supersonic 
flows by simply interchanging the sign of the parameter a. 


CALCULATION OF THE PARTICULAR INTEGRAL OF THE THIRD-ORDER 

ITERATION EQUATION 


In this section, the particular integral of equation (9) involving 
only <&2> their derivatives is derived. For this purpose, 

the variables x,y and the parameter a are introduced. Equation (9) 
then takes the following form: 


®3xx + $ 3yy = 


j(l + <O(*:ba02x +®1x®2xx) + [ 2 3 2 (! + <0 - l] $ lx $ ly $ lxy 
(1 + cr)|s0 2 (l + ff) - ^lxx^ix 2 + tK^P 2 “ l) $ lxx $ ly 2 + 


$ lxy $ 2y + $ ly t 2xy 


(19) 
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Use is again made of the complex conjugate variables z and z as 
independent variables. Thus, 


$2 = K.P.V2(z,z) 


$2x - R.P.(w 2z + v 2z ) 

®2y = R.P.i(v 2z - v 2 -) 

®2xx = R.P.(w 2 2 2 + 2w 2z5 + w 2 ~) 

®2xy = R*R*i( w 2zz ~ w 2zz^ 

® IXX^lX = R.P-lj:(«lz + ' ^15)^122 
*l*Ay 2 = - R - P ‘i (w lz - "lz^izz 

$ lxx $ 2x + #lx*2xx “ E - p -|( W lzz + W lzz)( w 2z + w 2z) + 

l 

R.P.ICw^ + wi^)(w 2zz + 2 w 2zz + v 2z £) 

, ®lxy*2y “ -R.P.§( v lzz, " hzz^ v 2z “ w 2z) 

$ ly®2xy = _ ■ R,P, 2^ w lz “ ^lz^ w 2zz ” w 2zz^ 
®Lx®ly^lxy = " R * P, ^( w lz - %z^) w lzz 
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Then equation ( 19 ), with $3 = S.P.w 3 (z,z), can he written as follows: 
w 3zz = " P 2 ) a ( w 2z w lzz + v 2zz w lz + v 2z w lzz + w 2zz w lz) + 

|(1 - P 2 )(2 + *)( w 2 z*lzz + W 2zz*lz + W 2z W lzz + *222^ + 

1(1 - p2) a [ p 2 (1 + 2 a) -,(1 - 2a)]w lz 2v lzz + 

Jg(l _ p2jjp2(4 + 5cr + 2a 2 ) + a (3 + 2a)J w^g^lzz + 

1(1 - P 2 )[p 2 (2 + 5a + 2a 2 ) + (-2 + a + 2a 2 )] v lz w lz w lzz 

introducing the expressions for Hfe(z.S) allowing: 

respect to z and 5 from equation ( 15 ) yields for w 3z g the folio g 

w 3z - = fel - P 2 ) 2 a 2 5(w lz 2 w lzz ) z + ^(1 - P 2 ) 2 *(2 + ^)(^l)z( w lz w lzz)z + 


3^(1 - p 2 ) 2 a(w lz F z ) z + ^(1 - P 2 )d[p 2 (2 + a) + (b + 3a] w lg (w lz 2 ) z + 


~(l - P 2 ) 2 a(2 + a)(w lz 2 ) z (zw lz ) z + 


^(1 - P 2 ) 2 (2 + °) 2 v lzz (vihzh + 


^(1 - p2)2(2 + *)(* 15 F z ) Z z + - » 2 > 2 < 2 + s ) 2 (0 *»Uk * 

i(l - p 2 )a 2 (5 + 3P 2 )(w 1z 3 ) z + 

96 

i(l - p2)[p2(8 + 10a + 3 a 2 ) + 0(6 + 5 a)] %zz w lz 2 


( 20 ) 
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In the derivation of this expression for w^g, free use was Bi&cLe 
of the fact that insofar as the real part $3 of W3 is concerned, 
terms of the nature g(z)£(z) and g(z)h(z) are equivalent. It is 
important to note that this type of operation leaves unaltered the real 
part °f V3. Since $3 is the quantity sought for in the calcu- 

lations, changes in the imaginary part of V3 are of no consequence 
in the final results. 

Equation (20) can he integrated immediately by inspection and yields 
the following result: 

W 3 = i (l " p2)2 ° 252w lz 2w lzz + ^( 1 - 3 2 ) 2 ct5w 1z f z + 

^(1 - P 2 ) 2 a(2 + ff)ztr 1 w lz w l22 + 

p 2 )a[p 2 (2 + a) + (4 + 3cr)J v-jW^ 2 + 
p 2 ) 2 a(2 + <j)zw l5 w lz 2 + jg(l - P 2 ) 2 (2 + tr) 2 w 1 ff 1 gw lz + 


3S (l " 

— (l - 
32' i 


Ig(l - P 2 ) 2 (2 + a)(Fi^ l5 + F z w x ) + ^(1 - p 2 ) 2 (2 + a) 2 w 1 2 w lzz + 


^(1 - P 2 )a 2 (3P 2 + 5)5 v 1z 3 + ^(1 - p 2 )[p 2 (8 + 10a + 3a 2 ) + 


a(6 + 5a 



dz 


( 21 ) 


Equation (21) is the particular integral of equation (20). The 
most general solution is obtained by adding an arbitrary analytic 
function G(z), satisfying the homogeneous or Laplace equation G zz = 0. 
An arbitrary function of z, customarily included in the general solution, 
need not be considered here because only the real part $3 of W3 is 

of interest. In fact, the omitted arbitrary function is the complex 
conjugate G(z). 
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In order to obtain tbe desired form of $3 ( ihe real part of W3) 

from equation (21), F(z) is replaced by its expression given in 
equation (l 6), and the real part of f(z) in equation (17) is replaced 

by 


R.P.f(z) = + 2oy^2x®ly - 2 (2 + crj^ix 

M 

00 

The final form of the particular integral of the third-order iteration 
equation (9) then becomes 


* 3 = - |(1 - 3 2 )oy(® lx f 2y + + l^ 1 - P 2 )( 2 + + *i* 2 x) + 

0(1 ‘ ~ ®ly 2 ) " " P 2 ) 2cj2 y 2 ®ix^y 4, Ixy + 


^(1 - p 2 )a[(6 + 5cr) + 3(-2- + 0)32] yd^y.3 + 

^(1 - 3 2 )ajj[lO - or) - (10 + Yajpsjy^. 2 ^ + 

I ( 1 ' p2)2ff(2 + a)y *l(*lx®Lxy + $ ly®Lxx) + 

J£(l - p 2 ) |(-l6 - lOcr + a 2 ) + (16 + 22a + Ttr 2 )? 2 ]^^ 2 
^(1 - P 2 )a |(6 + 5a) + 3(2 + ajp 2 ^^ 2 _ 

^(1 - p 2 ) (2 + cr) 2 ®^^ + J^(l - p 2 ) Qr(6 + 5a) + 

0 

(8 + 10a + 3^2)3^!^^ J'J(® lx 2 - $ ly 2 )dx + 2®^®^ dyj 


( 22 ) 
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The corresponding expression for for supersonic flow is 

* M 2 

7+1 oo 

obtained by simply replacing a by -a with a = —= and 

o Ho 2 - 1 

P^ by ~p2 with P^ = _ i. The physical, plane variables X and 

are easily inserted into both equations (18) and (22) by means of the 
transformation equations (10), 


x = X 
y = PY 


It is pointed out that the forms of the two particular integrals, 
equations (l8) and (22), derived in this paper are identical for both 
subsonic and' supersonic flow. The apparent differences are caused by 
a change in sign of the parameter p2. Thus, p2 and a are 
positive for both subsonic and supersonic flow. Actually, of course, 
the functions represented by <&^, . . . are different for 

the two types of flow. For subsonic flow, these functions are derived 
from analytic and nonanalytic functions of z and 2j whereas for 
supersonic flow, they involve the real "characteristics" variables x ± Py. 

Note that the last term of the expression on the right-hand side 
of equation (22) contains the indefinite integral 


r - + 24 lxV <*•] < 23 > 

It is obvious from the corresponding complex integral in equation (21) 
that the integrand of equation (23) is an exact differential. This fact 
can also be easily verified with the help of Laplace’s equation. 


^Lxx + ®lyy = 0 


Thus, 


l^ 2 - v) - y^v) 


(24) 
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Equation (24) represents the necessary and sufficient condition that the 
integrand of equation (23) "be an exact differential. Further, according 
to the theory of exact differential equations, the integral I may he 
expressed as follows: • 


where 



M =®1* 2 -»ly 2 
K = 2$ lx®ly 


(25) 


and where hy J' M dx is meant the result of integrating M dx with y 
considered constant. The expression within the parentheses, namely 


N 



is a function of y only. This statement can be verified as follows: 

Thus, 

i - h kT* 

and because y is considered constant in the process of integrating M dx, 
it is clear that 


U* M dx = M 


Hence the right-hand side of equation (26) is 
because of the condition for the existence of 


dN dM 
3x " 'Sy’ 


which vanishes 


an exact differential. 
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Note that In general it is simpler to perform the complex 


integration 


J v lz 2 


rather than transform to a areal 1 integral and 


then perform the integration. 


Langley Aeronautical'. Laboratory 

National Advisory Committee for Aeronautics 

Langley Air Force Base, Va., May 29, 1950 
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